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METHOD OP EXAMINING AN OBJECT BY MEANS OF ELASTOGRAPHY 



The invention relates to a method of examining an object by means of 
elastography, as well as to an arrangement for carrying out such a method and to a computer 
program for controlling such an arrangement. 

The MR elastography (MR = Magnetic Resonance) method which is disclosed 

5 in United States patent application 09/743659 (PHD 99-152) utilizes the fact that the phase in 
an MR image of the object changes under the influence of mechanical oscillations acting on 
the object. The extent of such a change is dependent on the excursion of the tissue under the 
influence of the mechanical oscillations. Information concerning given mechanical 
parameters of the tissue can thus be derived from MR phase images, that is, images 

10 reproducing the phase of the nuclear magnetization; for example, such a parameters relates to 
the elasticity. In ultrasound elastography even direct measurement of the excursion is 
possible. 

Malignant lesions such as, for example, a mastocarcinoma, clearly deviate in 
respect of elasticity from normal healthy tissue. However, there are also benign lesions which 
1 5 have similar elastic properties. Anomalous elastic properties, therefore, per se do not yet form 
an adequate indication that a tissue is malignant. 

It is an object of the present invention to provide a method of elastography 
such that additional information is derived which enables better differentiation between 
different types of tissue. This object is achieved by means of a method for examining an 
20 object in accordance with the invention which includes the steps of: 

a) exciting the object by means of mechanical oscillations, 

b) measuring the variation in time of the excursion from the rest state of voxels of the object 
which are subject to the waves caused by the oscillations in the object, 

c) determining the non-linear distortions from the variation in tune of the excursion, and 
25 d) evaluating the non-linear distortions. 

The mvention is based on the following considerations and insights. When a 
pressure exerted on biological tissue is varied, the deformation of the tissue does not vary to 
the same extent as the pressure, because the modulus of elasticity increases (the tissue 
becomes harder). This non-linear relationship between pressure and deformation is 
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particularly pronounced in the case of malignant lesions and is the reason why, when the 
tissue is excited by mechanical oscillations, the variation in time of the excursion no longer 
corresponds to the variation in time of the mechanical oscillations but exhibits (non-linear) 
distortions. Therefore, the non-linear distortions in the temporal variation of the excursion of 
a voxel are a measure of the non-linearity; such a measure is suitable for the classification of 
tissue. In the case of excitation by means of, for example, mechanical oscillations which vary 
sinusoidally in time, such distortions cause the excursion of the voxels in the object to be no 
longer purely sinusoidal and that it additionally exhibits higher harmonics of the mechanical 
oscillations, notably harmonics of twice the frequency. 

Each voxel could be individually measured so as to determine the non-linear 
distortions (and possibly also the elasticity) measured for each voxel. The version disclosed 
in claim 2, however, makes it easier for the examiner to form a diagnosis for a larger region. 

Claim 3 enables the invention to be carried out by means of MR elastography 
means, thus enabling location-dependent determination of the modulus of elasticity in 
addition to the non-linearity. In principle, however, the method can also be canied out by 
means of ultrasound. 

Claim 4 discloses a version of the method of claim 3 which enables the 
simultaneous measurement of excursions of a voxel at the fundamental frequency and at 
twice the fundamental frequency. Instead, however, it would also be possible to compose the 
magnetic gradient field from a component having the fundamental frequency and a 
component having twice the fundamental frequency. 

The further version as disclosed in claim 5 prevents field inhomogeneities 
from exerting a disturbing effect on the phases. In this case the two individual sinusoidal 
oscillations of the gradient field must have the same phase position. 

The version disclosed in claim 6 takes into account the fact that the (location- 
dependent) non-linearity and elasticity tissue enable improved classification of the tissue. The 
method disclosed in claim 3, however, does not require additional measurements for this 
purpose. 

Claim 7 defines an MR apparatus for carrying out the method in accordance 
with the invention and claim 8 discloses a computer program for controlling an MR apparatus 
of this kind. 

The invention will be described in detail hereinafter with reference to the 

drawings. Therein: 

Fig. 1 shows a diagram for the classification of different types of tissue. 
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Fig. 2 shows an MR apparatus for carrying out the method in accordance with 

the invention. 

Fig. 3 shows a flow chart illustrating the acquisition of MR phase images. 
Fig. 4 shows one of the sequences intended for the MR acquisition. 
Fig. 5 is a diagrammatic representation of a part of the evaluation process, and 
Fig. 6 shows a flow chart of the evaluation process. 

The following known relation exists between the pressure P on a sample and 
the deformation s of the sample due to this pressure: 

P = 8E 0) 

Therein, E is the elasticity of the tissue. For the deformation 8 it holds that: 

e = (L„-L)/Lo (2) 

Therein, LO is the extent of the non-compressed sample and L the extent of the compressed 
sample. As has already been stated, in the case of samples of a biological tissue the 
deformation 8 is not proportional to the pressure P, because the elasticity mcreases as the 
deformation increases. This increase can be suitably approximated by the equation 



(3) 



Therein, EO is the initial elasticity (that is, the elasticity without compression) and tiie 
exponent a is a measure of the non-linearity. For a = 0 a linear relationship exists between P 
and 8. The larger a, the more non-linear Ihis relationship will be. 

Fig. 1 shows a diagram in which the elasticity Eq is plotted as tiie ordmate and 
the exponent a, characterizing the non-linearity, is plotted as the abscissa. With each sample 
in Fig. 1 there can be associated a point which characterizes the elasticity Eo and the 
exponent a for this sample. In the diagram a number of tissue samples are marked by a 
respective symbol. The symbol "+" denotes tissue samples from a mastocarcinoma. The 
symbol "o" denotes benign lesions in tiie form of a fibroadenoma, the symbol "□" denotes fat 
tissue and the symbol "A" denotes tissue samples of growths in the milk ducts which could 
change from a benign to a malignant condition. 
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Tissue samples of mastocarcinomas have a higher elasticity (i.e. they are 
harder) and at the same time a more pronounced non-linear relationship between the pressure 
P and the deformation e in comparison with healthy tissue. Therefore, tissue samples can be 
classified on the basis of these parameters. Tissue samples situated to the right of the straight 
line in the diagram are indicative of a maUgnant lesion whereas tissue samples situated to the 
left thereof are indicative of a benign lesion. 

The diagram of Fig. 1 is based on in vitro examinations of tissue samples 
during which &e tissue samples were exposed to an increasing pressure P and the 
deformation s was measured. The examination of tissue of living patients in this manner is 
not without difficvdty. The mvention is based on the recognition of the fact that the same 
parameters can be non-invasively determined for living patients while utiUzing elastography 
means. When a living object is excited by stationary mechanical oscillations, the excursion of 
the voxels in the mechanical tissue should exhibit the same variation in time as the 
mechanical oscillations if the elasticity were constant and mdependent from the deformation. 
However, when a voxel is situated within a tissue which exhibits a non-linear relationship 
between the pressure P and the deformation s, the excursion of tiie voxel can no longer 
exactly follow the variation in time of the oscillations and distortions occur which are a 
measure of the non-linearity. The measurement of these distortions or non-linearity will be 
described m detail hereinafter on the basis of an example involving MR elastography. 

Fig. 2 is a diagrammatic sectional view of an MR apparatus for carrying out 
this method and also shows the hardware and software components required for the operation 
thereof. The MR apparatus includes a main field magnet 30 for generating a main magnetic 
field Bo which extends horizontally and parallel to tiie plane of drawing of Fig. 2. Inside tius 
cylindrical main field magnet there is arranged a set 3 1 of gradient coils which serve to 
generate magnetic fields having gradients in tiiree mutually perpendicular directions. Inside 
these gradient coils there is anranged an RF coil 32 for generating an RF magnetic field. 

The currents for the gradient coils 3 1 are supplied by a generator arrangement 
33 whereas the currents for the RF coil 32 are supplied by an RF generator 34. The variation 
in time of the currents generated by the generators 33 and 34 is determined and controlled by 
a control unit 35. The MR signals generated m the mammae during the examination are 
received by MR coils (not shown) so as to be further processed in a receiver 36. The output 
signals of the receiver 36 are ^plied to an image processing unit 37 which reconstructs an 
MR image from the MR signals received; this MR image is displayed on a monitor 38. 



wo 03/073933 PCT/IB03/00757 

5 

In the examination zone defined by the coils 3 1 and 32 there is arranged a 
female patient 47 who is positioned on a cushion 41 resting on a patient table 2. In firont of 
the cushion 41 there is arranged a mammography accessory 1 which has the same height as 
the cvishion 41. The mammae of the patient hang down in openings of the mammography 
accessory in which they are exposed to mechanical oscillations while being slightly 
compressed, said oscillations propagating in the horizontal direction (fi'om the feet to the 
head) while varying sinusoidally in time. The variation in time of the mechanical oscillations 
produced in the mammography accessory 1 are synchronized, by means of flie control unit 
35, with the variation in time of the cvirrents generated by the generators 33 and 34. 

The above MR apparatus is described in detail in German patent application 
10156178.4 (PHDE 010302) whereto reference is made m order to avoid duplication. 

Fig. 3 illustrates the execution in time of the acquisition of the MR data for an 
MR elastography examination by means of the MR apparatus shown in Fig. 2. After the 
completion of the mitialization (100) and after the patient 40 has been positioned m 
conformity with Fig. 1, the mammae of the patient are slightly compressed so that they are 
fixed and still elastically deformable. Subsequently, mechanical oscillations which vary 
sinusoidally in tune are generated in the mammography accessory 1 , said oscillations acting 
on the mammae of the patient and lasting for the entire MR acquisition (block 101). 

After said oscillations have become stationary, in block 102 there is first 
generated a first MR sequence whose variation in time is shown in Fig. 4. The sequence 
comprises a 90° RF pulse RFi and subsequently a refocusing 180* RF pulse RF2 (fourth line 
of Fig. 4). Between the two RF pulses RFi and RF2 there is generated a single, temporally 
sinusoidally varying oscillation Gi of the magnetic gradient field (third line of Fig. 4); the 
period duration of this oscillation amounts to exactly half (for example, 10 ms in the case of a 
fi-equency of the mechanical oscillation of 50 Hz) of the period duration of the mechanical 
oscillations. After the second RF pulse, a second oscillation G2 of the magnetic gradient field 
G is generated; this second oscillation has exactly the same variation in time as die oscillation 
Gi and its beginning has been shifted by one quarter period (5 ms) of the mechanical 
oscillation relative to the end of the first oscillation. 

During the two RF pulses RFi and RF2, a magnetic gradient field is applied in 
the z direction, said gradient field allowing the RF pulses to become active exclusively in a 
slice extending perpendicularly to the longitudinal du^ction of the patient; like the phase 
encoding gradient and the read gradient, this gradient field is not shown either in Fig. 4. 
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The first line of Fig. 4 shows a part of the variation in time of the mechanical 
oscillations and the second line shows a higher harmonic of these oscillations at twice the 
frequency. Such higher harmonics can occur continuously in regions in which a non-linear 
relationship exists between the deformation and the elasticity when the waves excited by the 
mechanical oscillations propagate in the mammae and reach such regions. Moreover, the 
dashed lines of these two first lines indicate that the spins deflected at the relevant firequency 
are refocused by the pulse RF2. 

The gradient oscillations Gi and G2 change the phase of the nuclear spins 
deflected by the mechanical oscillations. The extent of the change is dependent on the degree 
of excursion in the direction of the relevant gradient, whereas the sign of the change is 
dependent on the phase position between the gradient oscillation and the excursion of the 
nuclear spins. 

It appears firom fig. 4 that the sign of the gradient oscillation Gi changes 
whereas the sign of the mechanical 50 Hz oscillation (first line) remains the same. Therefore, 
no change of phase occurs. However, the second gradient oscillation G2 causes a change of 
phase because the 50 Hz oscillation has its zero crossing at the same instant as the gradient 
oscillation G2. The two gradient oscillations Gi and G2 cause equally large changes of phase 
for the spins deflected at the harmonic (100 Hz) of the fundamental firequency, said changes 
being summed because of the variations shown. 

The sequence shown in Fig. 4 is sensitive to the excursion at the fundamental 
firequency as well as to the excursion at twice the fundamental firequency. If this were not tiie 
case, the fimdamental fi«quencies of the excursion and their harmonic would have to be 
separately measured, so that the measuring time would be prolonged. 

After the reception of the resultant MR signal in the step 102, the MR 
sequence is repeated in the step 103 for a different sUce which has been shifl:ed in the z 
direction (extending horizontally in Fig. 2) (block 103). Merely the carrier frequency of the 
RF pulses RFi and RF2 is then varied in such a manner that the other slice is excited while 
the variations in time of the magnetic gradient fields remain the same. The loop comprising 
flie steps 102 and 103 is tibien repeated as many times as there are slices in the three- 
dimensional region to be examined, for example, 20 times, each time a different slice being 
excited. 

Subsequently, the phase encoding gradient is changed in the step 104 and all 
slices are excited again so as to acquire the spin echo signals generated therein. Each 
individual slice is then excited with the same position in tune of the RF pulses and the 
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magnetic gradient fields in relation to the mechanical oscillation. The phase encoding 
gradient is changed as many times as there are phase encoding steps required, for example, 
128 or 256 times. The MR data are thus acquured from all slices so that an MR image can be 
reconstructed therefrom; this image represents the phase of the nuclear magnetization 
distribution in the three-dimensional region formed by the slices. 

After the acquisition of the MR data for the first MR phase image of the three- 
dimensional region, m the step 105 a shift m time takes place between the MR sequences and 
the mechanical oscillations, that is, a shift such that the beginning of an MR experiment is 
shifted by a fraction of a period AT; for example, AT = T/8 (step 105). Subsequently, the 
loop comprising the steps 102, 103 and 104 is completed again so that a fiirther MR phase 
image of the three-dimensional region being examined is obtained. The two MR phase 
images differ only m respect of the regions in which the tissue or the nuclear spins excited 
therein have been deflected in the x direction. After that there are formed further phase 
images for which the shift in time between the mechanical oscillations and the sequences is 
changed (for example, AT = 2T/8,3T/8 ... 7T/8. The data entered by the user or fetched by the 
entry made by the user can be displayed on the monitor 8 or also on a separate display. 

Finally the MR data of a number of MR phase images (8 in the present 
example) have been acquired, said MR data being dependent on the excursion of the spins in 
the X direction, tihat is, assuming fliat the gradient of the gradient fields Gi, G2 extends in the 
X direction. 

In order to enable measurement of the excursion also in a direction other than 
the x direction, the sinusoidal gradients Gi and G2 are formed in a direction other than the x 
direction (block 106), for example, in the y direction, so that, after completion of the loop 
comprising the steps 102, 103, 104 and 105, there is obtained a set of MR phase images of 
the three-dimensional region which has been influenced by the excursion of the spins ui the y 
direction. Finally, this procedure is repeated for the z direction (that is, the sinusoidal 
gradients Gi and G2 are applied in the z direction), so that finally there are available three sets 
of 8 MR phase images each which are dependent on tiie excursion of the spins in the x 
direction, the y direction or the z direction. The acquisition method is then terminated (block 
107). 

The processing of the acquired MR data will be described in detail hereinafter 
witifci reference to Fig. 6. As has akeady been stated, after the initialization in fbe block 200 
MR phase images are formed from the acquired MR data in the step 201 . In the step 202 the 
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variation in time of the excursion Ux, Uy and Uz is determined for each of tibie three spatial 
directions, that is, for the individual voxels of the three-dimensional examination zone. 

For the pxirpose of illustration reference is made to the diagrammatic 
representation of Fig. 5. The upper line diagrammatically shows the three sets of MR phase 
images which are dependent on the excursion in the x direction, the y dkection and the z 
direction, respectively. In each MR image associated with a set the phase is considered in the 
same voxel. The eight images yield eight base points for the variation in time of the phase or 
the excursion in the x durection, the y direction or the z direction in the relevant voxel. When 
this voxel is situated in a region exhibiting non-linear elastic properties, the excursion in this 
voxel will not be purely sinusoidal in time, because the deflection in this voxel contains, in 
addition to a component having the frequency of the fundamental oscillation, also higher 
harmonics, notably a component having twice the fundamental frequency. 

In order to determine the amplitude of these components, in the step 203 a 
Fourier transformation is carried out over the phase or the excursion of the voxel as defined 
by the eight base points. The Fovirier transformation yields three quantities aix, aiy and aiz 
which represent the amplitudes of the excursion in the x direction, the y direction and the z 
direction at the fundamental frequency, and also yields three further quantities a2x5 aiy and aiz 
which represent the amplitudes of the excvirsion m the x direction, the y dkection and the z 
direction at double the fundamental frequency. The amplitude ai of the component having the 
fundamental frequency can be calculated therefrom in conformity with the equation: 



a,=.j4+4+4 (4) 

For the amplitude of the component having twice the fundamental frequency it 



holds tiiat: 



When the two components are related to one another in conformity with the 



equation 



(6) 
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there will be obtained a value a which is a measure of the non-linearity. Granted, the 
amplitude 32 per se would already be a measure of the non-linearity, but a2 is dependent not 
only on the non-linearity but also on the attenuation of the wave produced by the excitation in 

5 the object. The parameter a is independent therefrom to a high degree. It corresponds to the 
value a of the non-linearity in conformity with the equation 3. This correspondence is given 
via the exact application of the distortion tensor as well as the free energy while taking into 
account higher terms. The relationship between a and a increases strictly monotonously. 
Therefore, analogously to Fig. 1 , an equivalent diagnostic opinion can be given by way of Eo 

10 and a. 

After a respective value a has been determined for all voxels, in the step 204 
images can be formed which represent the non-linearity of the elasticity in the examination 
zone. 

Because the non-linearity per se generally does not allow an opinion to be 

1 5 expressed as regards the benign or malignant nature of a lesion, as described with reference 
to Fig. 1, in the step 205 there is furthermore generated an image which represents the spatial 
distribution of the modulus of elasticity in the examination zone. The measuring values 
required for tiiis purpose are available ahready after the acquisition in conformity with Fig. 3. 
The calculation of the modulus of elasticity from these values is described in detail in the 

20 previously mentioned United States patent application 09/743659 (PHD 99. 1 52) which is 
hereby specifically referred to. 

The user can then select individual regions in an image, for example, regions 
marked by means of a cursor, and be presented with the values of E and a for these regions, 
for example, in a diagram analogous to Fig. 1 (step 206). The method is then terminated. 

25 As has already been stated, ultrasoxmd elastography means are also suitable for 

tracking the propagation of the waves excited by the mechanical oscillations in the object, 
such tracking being significantly faster than in MR elastography. The non-linear distortion of 
the variation in time, resulting from the non-linear properties of the tissue, can then be 
determined in a manner analogous to that described above. In addition to the mammae, other 

30 regions may also be examined, for example, the prostate. 



